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ABSTRACT 



The algorithm presented is an extension of the Land and 
Dolg branch and bound method combined with the branch 
selection techniques presented by Beale and Small to solve 
integer or mi xed- I nteger linear programs. The algorithm 
obtains the solution by solving a linear program with upper 
and/or lower bounds on selected branch variables. By 
systematically changing these bounds# and maintaining only 
the current canonical form# the solution is assured using a 
minimum of excess computer storage above that required to 
solve the linear programming problem. Thus the problem can 
be solved entirely within the computer core# and the problem 
converges to the solution faster than most other general 
integer linear programming algorithms. 
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1. INTRODUCTION 



In the literature today, there are many algorithms to 
solve Integer or ml xed- Integer linear programs, but most of 
these algorithms are for specialized types of problems, such 
as transportation problems or machine scheduling problems. 
As Bellmore and Nemhauser CO say, "Integer linear 
programming algorithms are notorious for converging rapidly 
on one problem and then performing miserably on the next." 

The algorithm presented here Is Harold Greenberg's 
natural extension of the Land and Doig branch and bound 
method CO , E| combined with branch selection techniques 
presented by E.M.L. Beale and R.E. Small CO , to solve 
Integer Linear Programs. In addition, the algorithm has 
been developed to solve mi xed- i nteger linear programs where 
some, but not all of the variables are i nteger- rest r i cted. 

The problem is solved by primal linear programming 
using upper and lower bounds on an intelligently selected 
integer-restricted variable that does not have an integer 
value in the current solution. These bounds are system- 
atically changed until an optimal solution is assured. In 
addition to speeding convergence by having less branches 
than most other methods, this method also uses a minimum of 
excess computer storage above that required to solve the 
linear programming problem. Thus the need for offline 
storage is eliminated and the entire problem is solved and 
maintained in the core, giving a relatively small cpu time. 
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Computational results and 
are included to show this 



compar i sons 
small time. 



with other algorithms 
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LL FORMULATION OF THE PROBLEM 



The problem to be solved Is: find X * (x^, 

minimizes CX' 
when AX'=B 



,x n ) that 
(1) 



0 1 xj £ uj , j*l, . . . ,n 

Xj integer for some or all j*l,...,n. 

C is an n-vector, B is an m-vector, and A is an mxn matrix, 
and the Uj are positive integers (or are infinite). 

To formulate a problem into this form, one need not 
restrict the slack/surplus variables to being integer, thus 
there is no need to restrict the elements of the A, B, or C 
vectors to being integer. 

The solution to (1) includes the case where Xj is 
restricted to being 0 or 1. This accomplished by 

restricting Xj to being integer with an upper bound of one 
(i.e. the corresponding u. equals one). 
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III. THE BRANCHING PROCEDURE 



A. THE LAND AND DOIG METHOD 

The branching procedure employed In thls> algorithm Is 
an extension of the branch and bound method of Land and Doig 
|9|; their method will be explained first. Land and Doig 
first solve (1) as a linear program to obtain the continuous 
solution, (z° , x°,...,x°). If the continuous solution Is 
all Integer, the problem Is solved. If there are fractions, 
this Initial value of the objective function, z°, serves as 
a lower bound on the value of any feasible optimal 
Integer-solution. One variable, Xj, Is then selected whose 
value Is not Integer. In the final integer-solution, xj , 
must satlsify either Xj <. [x° ], or xj J> [x° ] + 1, where 
Cx° ] Is the greatest Integer less than Xj . The procedure 
initially Involves the solution of two problems: both 
contain the Initial linear programs, but with an additional 
constraint; In one case the constraint Is Xj * [xj ], and 
In the other case It Is Xj * C*j 3+1. Once these 
solutions are obtained, and supposing that both are not 
Integer-solutions, the one with the minimum objective value 
Is investigated first, while the other Is stored for later 
Investigation. Suppose the minimum of the two objective 
functions Is z* , with solution: 
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whose 



(z 1 ,x | , . . . , [x ?J , . . . ,x 1 ) . Now a new variable, x. , 
i j n k 

value Is not Integer, Is selected, and three new problems 
are solved: these are the Initial linear program (1) subject 
to the additional constraints Xj * [xj ], x^ = [x£ }, In the 
first case, Xj ■ [x° J, x^ » J ♦ 1/ In the second 

case, and Xj = [xj 1 ] - 1, In the third case, all assuming 
CXjl^O. 

The procedure Is again repeated until an Integer- 

solution Is obtained. At every step where there Is a new 

variable chosen for restriction, there Is an additional 

problem to solve. Thus the next step would have four 

problems to solve, and the one after that, five. If the 

Integer-solution Is not obtained. Once an integer-solution 

q 

Is achieved. Its objective function value, z , becomes a 
new upper bound on the optimal solution. With this bound, 
the procedure then Investigates the other problems that are 
solved at this step. At any point in the branching, when 
the current upper bound is exceeded, or when an 
Integer-solution Is obtained, that branch is terminated. 

This complete process enumerates the possibilities for 
the solution to (1) In a directed tree. The rooted node 
corresponds to z° and It has directed branches to the nodes 
corresponding to the solution of the linear program (1), 
with Xj set at a specific value, each of these Integer 
values representing a single branch. At each of these 
nodes, the branching Is repeated. The minimum objective 
value at each node determines the next branch, which limits 
the potentially large size of the tree. 
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The problems that are stored for further Investigation 
rapidly create an immense storage problem/ thus the Land and 
Doig method becomes unwieldy for problems with more than 
Just a few variables. 

The algorithm presented here takes the Initial thought 
of the Land and Doig method, but uses an Inequality 
constraint for the selected variable. This Inequality may 
be as an upper or a lower bound for the variable. By 
employing a bounded variable linear programming solution 
technique similar to that from Dantzig [ 3 ], the algorithm 
allows the problem to be solved with a minimum of excess 
storage used above that storage required to solve the linear 
programming problem. 



B. THE NEW BRANCHING PROCEDURE 

As with the Land and Doig method, the procedure obtains 
the continuous solution to the problem. If this solution is 
not an integer solution, the procedure then selects a 
non-integer valued variable, Xj , to branch upon. This 
creates two new problems; both contain the original linear 
program Cl), but one has the added constraint of Xj<.Qx^], 
while the other has the added constraint of Xj >. Cx?] + 1. 
Initially, only one of these problems is solved. The 
solution to this problem is treated in the same manner, and 
the process Is continued until an Integer solution is 
obtained, or the problem becomes Infeasible. One thing to 
notice is at every node there are only two possible 
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branches, and there are no problems saved. This limiting of 
branches, alone, tends to reduce the number of nodes 
required for the solution of the problem. In addition, once 
an Integer solution Is obtained. It serves as an upper bound 
on the final solution, and no branching continues beyond a 
node where the objective value excedes the current upper 
bound solution. 

The strategy used Is to Investigate a single branch of 
the tree until z q >. zz , where zz Is the current upper 
bound on the optimal solution. At this point, backtracking 
begins, to find a node that has only one previous branching. 
This backtracking Is accomplished simply by removing all 
bounds from the nodes that are further out on the tree, and 
solving the resulting problem. This solution Is the 
re-creation of a previous solution that would have been 
stored otherwise. The second branch Is then made utilizing 
the Initial problem at that node which Is re-created from 
the latest transformation of the equations, and the 
constraint which was Initially Ignored at that node, and 
solving this linear program. At times, the procedure may 
Impose further bounds on a previously bounded variable. 
This presents no problem since the Initial bounds will be 
sat Is l f led. 

Thus, this procedure stores only the current solution 
and the current bounds. The optimal solution Is the 
existing upper bound solution when It becomes Impossible to 
branch any further. This Is realized when backtracking to 
the rooted node a second time during the solution procedure. 
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IV. SIMPLEX PROCEDURES USED 



The program utilized the two phase simplex method to 
solve (1), together with a bounded variable technique 
similar to that In [33. By using this technique, the 

problem maintains m constraint equations throughout the 
solution procedure and has a solution where some of the 
variables will be at their upper bounds, some of the 
variables will be basic, and some of the variables will be 
at their lower bounds. 

At each node In the procedure, we are Interested In the 

solution to (1) which satlslfles the new bounds given by 

0 <. 1 . i. x. <. u. <. u., j=l,...,n, where u Is the problem 

J J J J j 

upper bound, which may be Infinite, Uj Is the Imposed upper 
bound generated by the branching procedure, and lj Is the 
Imposed lower bound, also generated by the branching 
procedure. After achieving a solution for the x^ which 
satlslfles the constraints and the Imposed bounds, we can 
obtain the following canonical form [3], from the constraint 
equations In (1): 



x 



1 



- £ 

j =m+l 



y • . x . 
y U J 



-z 



n 

♦ t 

j *m+l 



x. 

J 



d • , l -1, . . . ,m. 



(3) 



o 



-z / 



where the x. , 1=1, ...,m, are the basic variables, and the 
x , J*m+l,...,n, are the non-baslc variables. When there Is 

j 

a feasible solution to this problem (3), with upper and 
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lower hounds, we can use a variation of the usual simplex 
method to obtain optimality. This variation Is suggested hy 
utilizing a method to Improve a feasible solution that is 
not optimal. To Improve this feasible solution, we can 
Increase a non-baslc variable x at Its lower bound when 

j 

C < 0, or we can decrease a non-baslc variable x at I ts 

J j 

upper bound when c. >0, where c. Is the reduced cost 

J J 

coefficient for the j ^ variable as found In the objective 
function In the cononlcal form. These conditions for 

optimality come from the following theorem |5|: 

Theorem: 

any values of x , j=l,...,n, satislfylng the constraint 

j 

equations In (3) and 0 <. lj <. Xj <. Uj <. Uj, j=l,...,n, 
are an optimal solution for minimum z, with objective 
value z , Iff c >.0 for variables at their lower bound, 

°_ j 

1 , and c <.0 for variables at their upper bound, u . 

J j J 

The proof Is obvious, for any Increase In a variable at Its 



1 owe r 


bound, or 


any decrease 


In a va r 1 ab 1 e 


at 1 ts upper 


bound 


an only 


Increase z. 


More deta I 1 ed 


proofs may be 


found 


l rt 1 3 | , or 


in 16|. 
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V. THE VARIABLE SELECTION PROCEDURE 



While any non-integer variable In the present solution 
can be chosen for the branching, the size of the directed 
tree can be further reduced by selecting the one which will 
produce the optimal solution quickest, thus reducing the 
number of branchings and nodes. The method employed In this 
algorithm was developed by Beale and Small [4] , using a 
suggestion of Drlebeck. 

The process Involves comparisons of the marginal cost of 
Increasing and of decreasing each of the non-integer basic 
variables In the current solution, as In the cost ranging 
procedure of linear programming, where the marginal cost of 
changing each of the non-baslc variables Is considered, and 
the non-baslc variable with the greatest change which still 
allow feasibility Is entered Into the basis. The marginal 
cost,J)^., of Increasing x. is given by 



J> = Min 

Ul j,y < 0 J 
' j 



(C. / I y - j I ) , 



and the marginal cost,!)^., of decreasing x. Is given by 

D i i 



J> n . - Min 

Dl j / y j j > o J 



(c. / y. . ) . 
•J 



th 



where "cj , Is the reduced cost coeflclent for the j 
variable as found In the objective function In the canonical 
form. By applying the strategy of taking the variable that 
makes the worst of the two alternatives as bad as possible, 
we have the alternative that Is the most unlikely to contain 
the final solution, and we branch on this variable In the 
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other direction 



This strategy Is Impl Imented by computing^, . xf. 

L) i D i 1 

and JL y j = J)yj x d ** f. ), for each non-integer basic 
variable that Is restricted to being Integer. Where f. Is 

i 

the fractional part of x. , x. * [x. ] + f. . The branching Is 
then made on the variable that makes the larger of these two 
quantities as large as possible. Say that 

M f x -L D | > 

Max,/ . , we then go back to find I such that^ = Max J 
i Ul D i | D i 

Thus, the 1 ^ basic variable would be the one branched upon, 
and the direction of the branching Is x. > [x. 1+1, since 

i i 

It Is the^ and not that Is the greater. 
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VI. THE ALGORITHM 



A. THE BRANCHING PROCEDURE 

The branching procedure Is contained In the following 
algorithm: 

1. Def I nl tl ons 

S(zz,x, , . . . ,x ) Is the current upper bound solution 
1 n 

which represents a feasible solution to (1) with objective 
value zz, and having all variables at Integer values which 
are required to be Integer In the solution. 

I ( L ) Is the Index of the L** 1 bounded variable. 

t h 

B(L) Is the bound restricting the L bounded 
va r I abl e . 

a . L 

BSIGN(L) Indicates whether the bound on the L 
bounded variable Is an upper or lower bound: BSIGN(L) = + 1 

for an upper bound, and BSIGN(L) = - 1, for a lower bound. 

N(L) Is the counter to show the number of branchings 
for the L t * 1 bounded variable; Initially, N(L) = 0. 

2_, Logic Procedure 

a. Set L=0, and take zz as Infinite, unless there 
Is some known feasible Integer-solution to (1), to get the 
Initial S(zz,x 1 ,...,x n >. Solve (1) as a linear program to 
obtain the continuous solution. If this solution has 
Integer values for all of those variables that are Integer 
restricted, the problem Is solved. Otherwise go to b. 
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maintaining the canonical form of the solution to (1). 
b. Set L*L ♦ 1. Go to b(l). 

(1) Select the basic variable, x® , that is 

J 

fractional, but required to be Integer, which is the most 
unlikely to produce the optimal solution, and go to b(2). 

(2) Set BSIGN(L) opposite the way which will 
produce the most unlikely chance for the optimal solution 
and go to b(3) or b(4). 

(3) If BSIGN(L) - ♦ 1, set B ( L) = x° . Go to 

J 

c. 

(4) If BSIGN(L) - - 1, Set B ( L ) = x° + 1. Go 

to c. 



c. Solve the linear program using the maintained 
canonical form with the new bound on the variable x^ , 
where J = l(L). Go to the following cases, one of which 
will hold: 



(1) If the solution produces an objective 

value, z, with z > zz, go to d. 

(2) If the solution produces an objective value 

z, with z < zz, and the solution has Integers for all of the 

required variables, redefine S(zz,x, , . . . ,x ) as a new 

l n 

feasible Integer solution upper bound, where zz * z, x^ 

,...,x , is the new solution. Go to d. 
n 

(3) If the solution produces an objective 

value z, with z < zz, and the solution does not have 
integers for all of the integer-restricted variables, go to 
b. 
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(4) If the problem has an Infeasible solution. 



go to d. 

d. Set LL = L, and to to d(l). 

(1) If N( L) =2, go to d( 2 ) . Otherwise N(L)=1; 

go to d ( 3 ) . 

(2) If L=l, the current feasible solution 

S(zz,x ,...,x )is optimal, and the algorithm stops. 

1 n 

Otherwise, set L=L - 1 and go to d(l). 

(3) Solve the linear programming problem 

starting from the current canonical form with the bounds 
B(L), B(L + 1 ),..., B( LL) , removed. If j = l(L), and Xj does 
not become a basic variable, the solution Is non-unique; Xj 
Is then made a basic variable. Go to d(4) or d(5). 

(4) If BSIGN(L) = +1, set B(L) = B(L) +1, 

BSIGN(L) = -1, and N(L) = 2, and go to c. 

(5) If BSIGN(L) = -1, set B ( L ) = B(L) -1, 

BSIGN(L) = +1, and N(L) * 2, then go to c. 

This completes the branching part of the algorithm. 

B. THE VARIABLE SELECTION PROCEDURE 

The selection procedure for which basic variable to 
branch upon, and which way to branch Is contained In the 
following algorithm. 

1 . Def l n l t l ons 

I Is the l basic variable. 

J Is the j t ^ 1 column In the basis. 

I B ( I ) Is which problem variable Is the l*"^ basic 



va r i abl e 
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I MU Is the Index for which I has the max^^ • 

IML Is the Index for which I has the . 

DU(J) Is the current ratio * |Cj/yjj|,lf y j j < 0. 

D L ( J ) Is the current ratio = | c j / y j j I , I f y | j > 0. 

XDU Is the current value of the min DU(J). 

J 

XDL Is the current value of the mjn DL(J). 

XLDU Is the current value of the maxi. 

i Ul 

XLDL Is the current value of the 

BSIGN Is the Indicator for whether the branch 
requires an upper bound or a lower bound on the variable to 
be restricted. BSIGN = +1, for an upper bound, and BSIGN = 
-1, for a lower bound. 

M Is the number of constraint equations, I.e. the 
size of the basis. 

2. Logic Procedure 



variable 
or I f x 
XDU=XDL= 



a. Set XLDU=XLDL=IMU=IML= 0, 1= 0, and go to b. 

b. Set 1=1 + 1. If I B ( I ) Is an artificial 

, or If x |g(|) not restricted to being Integer, 
I B C I ) I S t nte ser, go back to a. Otherwise, set J = 0, 

Infinity, and go to c. 



c. Set J = J + 1. If c. = 0, repeat step c. If 

J 

y.j< 0, go to c(l). If y.^ > 0, go to c(2), otherwise, y^ = 
0, and repeat step c. 

(1) DU( J) = Icj/y.jl. If DU ( J ) < XDU, set XDU 
= DU(J), and go to c(3). Otherwise, go straight to c(3), 
without setting these values. 



C 2 ) DL ( J ) » | c j/y j j | . If DL( J ) < XDL, set XDL 
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a DL(J), and go to c(3). Otherwise, go straight to c(3), 
without setting these values. 



(3) If J < M, go to c, otherwise this Is the 

last comparison for the I basic variable and go to d. 

d. Go to the following cases. 

(1) If XDL Is still equal to Infinity, XLDL Is 
unchanged, otherwise. If XDL x fjg(j) > XLDL, set XLDL = 

XDL x f| and set IML = I, go to d(2). 

(2) If XDU is still equal to Infinity, XLDU Is 
unchanged, otherwise. If XDU x (1 - f^^ ) > XLDU, set XLDU 
= XDU x (1 -f , .), and set IMU = I, go to d(3). 

(3) If I > M, go to b, otherwise this Is the 

last basic variable to consider, and go to e. 

e. One of the following cases will hold. 

(1) If IMU and IML both equal zero, this 

procedure was called In error, and there are no variables to 

branch upon. 



(2) If XLDU > XLDL, then the variable to branch 

upon is IMU, and BSIGN = + 1, to set an upper bound on x . 

I MU 

(3) If XLDU < XLDL, then the variable to branch 

upon Is IML, and BSIGN = - 1, to set an lower bound on x,.., . 

I ML 

This completes the variable selection portion of the 



algorithm. 
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ILL, EXA M PLE 



The basic branching method, along with the selection 
procedure for which variable to branch upon are illustrated 
in the following mi xed- i nteger linear programming problem. 
This is the same problem used in Cl] to Illustrate the Land 
and Dolg procedure. 

Min 4x^ + 5x2 

when 3x^ + X 2 - x^ -2 

X 1 + ** x 2 “ x 4 =5 ( 4 ) 

3x 1 ♦ 2x 2 -x 5 = 7 

Xj >. 0, and integer, 1 * 1,...,5. 

The continuous solution produces the canonical form: 



x 



1 





= 12 . 

10 






i x + l x 

1Q X ^ 10 5 10 



x 3 





42 

10 



(5) 



-2 



+ 



l'0 X 4 




il2 
10 * 



Since the solution has fractions for at least one of the 
integer restricted variables, we use the selection procedure 
to tell which variable to branch upon, and which way to 
branch. 
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The computations to select the branching variable, and 
the direction to branch are Included for this first 
branch l ng : 

1=1 f = 0.8, 1 - f = 0.2. 

J = 1 

c" = 0, thus no calculation. 

j 

J = 2 



DL ( 2 ) 


= 

0.2 


= 3.5 


XDL 


II 

• 

LTI 

• 




DU ( 3 ) 


1.1 

0.4 


= 2.75 


XDU 


= 2.75. 





XLDL = 3.5 x 0.8 = 2.8; I ML = 1. 

XLDU = 2.75 x 0.2 = 0.55; IMU = 1. 

1=2 f = 0.8, 1 - f = 0.2. 

J = 1 

c. = 0, thus no calculations. 

J = 2 

DU(2) = = 2.33, 

0.3 

XDU = 2.33. 

J = 3 

DL ( 2 ) = -i-*- 1 = 11.0, 

0.1 

XDL = 11.0. 

XDL x f = 11.0 x 0.8 = 8.8 > XLDL. 

XDU x (1-f) = 2.33 x 0.2 = 0.467 < XLDU. 
XLDL = 8.8; I ML = 2. 

XLDU = 0.55; IMU = 1. 
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I - 3 



f = 0.2, 1 - f = 0.8 



J = 1 

c. = 0, thus no calculations. 

J = 2 

DL ( 3 ) = — = 2.33, 

0.3 

XDL = 2.33. 

J = 3 

DUO) = — = 1.0, 

X • X 

XDU = 1.0. 

XDL x f = 2.33 x 0.2 = 0.467 < XLDL . 

XDU x (1-f) = 1.0 x 0.8 = 0.8 > XLDU . 

XLDL = 8.8; IML = 2. 

XLDU = 0.8; IMU * 3. 

XLDL = 8.8 > 0.8 = XLDU, 
therefore: BSIGN = - 1, 

IM = IML * 2. 

Thus we Impose the lower bound on x^>. 1, and solve the 
problem (5) with this added constraint to get the canonical 
form for node 1: 



x + -x 
1 3 2 



1 

-x 
3 5 



7 

3 



x. + 2x „ 
4. 2 



1 

; x 5 



_8 

1 



X 3 X 2 X 5 



z ♦ 3 X 2 + 3 X 5 



28 

y 



( 6 ) 
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with (z,Xj, x 2 / Xj , , Xg ) * (35/3,5/3,1,4,2/3,0). 

Since some of the integer-restricted variables are still 
fractional, we again use the selection procedure to get the 
new bound x^ 2 2 , and solve the new problem taken from ( 6 ) 
with the added constraint, x^ 22 , to get the canonical form 
for node 2 : 




x, - x 3x, » - 2 (7) 

3 2 1 

-z + 5 x 2 + 4Xj ■ 0, 

where (z, x^, x 2 , x x^, x,. ) = (13,2,1,5,1,1), which is integer 
for all of the integer restricted variables, thus, this 
solution becomes the solution vector, S( zz, x^, . . . , x R ) . Next 
the bound x^ 2 2 , is removed, and the simplex procedure is 
applied to backtrack to the previous node to obtain ( 6 ) 
again. Now we impose the bound Xj 2 1, and solve ( 6 ) with 
this new constraint to get the canonical form for node 3 : 

x 4 + 5x i - 2x 5 = 9 



x + 
2 



K - x -x 

2 1 2 5 




x 

3 






( 8 ) 



-z 




— 14, 



where (z,x ,x 2 ,x , x ,x ) 



(14,1,2,3,4,0). 



Since z > 13, we 



again backtrack to a node where there is only one branch, in 
this case to node 0 , removing the bounds x^i. 1 , and x ^2 1 / 
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and solving by the simplex procedure to again get the 



canonical form (5). Then a new constraint of x^ £ 0, is 

added to (6) and the resulting problem is solved to get the 

canonical form for node 4: 

x ♦ 10x„ - 3x, = 8 

5 2 4 

x ♦ 4x - x * 5 

12 4 

x +llx - x « 13 (9) 



3 2 4 

-z -llx + 4x *- 20/ 

2 4 

where (z,x ,x ,x ,x ,x ) * (20/5/0/13/0/8), which again 

1 2 3 4 5 

has z > 13. Since there are no nodes with just one branch 

which can be backtracked to, the problem is solved and the 
solution is: (ZZ/X^/X^ ,x^,x^ ,x^) a (15,2/1/5/1/1). The 

procedure is exhibited in the tree shown in figure 1. 

This procedure needed only two nodes to obtain the 
solution, and four nodes to assure the solution, while this 
problem (4), as shown in [2] requires five nodes. Thus the 
advantages of this procedure over the Land and Doig 

procedure is evident even in this small problem. 
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Node 0 



z » 11.2 




Node 4 



Node 1 



z = 20 
z > zz 



z = 11.67 



x 



Node 3 



1 ' 




Node 2 



> 2 



z = 14 
z > zz 



z = 13 



( I nteger Sol ) 
zz = z = 13 



Node 0: (112/10,18/10,8/10,42/10,0,0) 

Node Is (35/3,5/3,1,4,2/3,0) 

Node 2: (13,2,1,5,1,1) 

Node 3: (14,1,2,3,4,0) 

Node 4: (20,5,0,13,0,8) 



Had we just used some arbitrary rule for selecting which 
non-integer basic variable to branch upon, the number of 
nodes may have been much greater. The following tree. In 
Figure 2, Illustrates the arbitrary rule of selecting the 
first non-integer variable In the basis, x^, and then 
branching above this variable by setting the bound. 



Figure 1 
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Figure 2 
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VIII. COMPUTATIONAL RESULTS 



The algorithm has been programmed in Fortran IV, and run 
on an IBM 360/C7. The complete program is given in the 
appendix. The computer core needed to run a problem is 62K 
plus 36n + 32m + 4m x (m + n) bytes, where n is the number 
of variables, and m is the number of constraints in the 
particular problem when in the form as (1). The program 
utilizes a simplex procedure in real variables with a 
round-off correction which requires that all problems using 

this program have all integer coefficients when in the form 
as ( 1 ) . 

The results in Figure 3 are for problems in [7], and the 
data for the other programs are from [7] and [8] . 

The results show that this branch and bound algorithm is 
faster than most of the general algorithms for integer or 
mi xed- i nteger linear programs. 
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Comparison of Computational Experience on Haldi and IBM 



Source 


m 

(g) 


n 

(h) 


LI PI 
(a) 


1 PM 2 
(b) 


1 PM 3 
(c) 


BBA 

(d) 


BVA 

(e) 


1 ter 
(f ) 


1 ter 


1 ter 


1 ter 


Time 
( sec ) 


T ime( j ) 
( sec ) 


Haldi 


















1 


4 


5 


25 


19 


51 


94 


3.15 


— 


2 


4 


5 


21 


18 


58 


26 


.83 


-- 


3 


4 


5 


22 


17 


80 


41 


1.33 


— 


4 


4 


5 


17 


29 


35 


12 


.38 


— 


5 


6 


5 


181 


765 


F(i ) 


22 




30 


6 


6 


5 


126 


20 


F 


245 


8.2 


12 


7 


4 


5 


173 


731 


F 


31 


1.03 


24 


8 


4 


5 


125 


20 


F 


31 


1.02 


9 


9 


6 


6 


43 


39 


192 


22 




8 


10 


10 


12 


120 


F 


F 






17 


11 


21 


56 


104 


— 


-- 






-- 


12 


21 


56 


103 


-- 


— 






— 


13 


21 


56 


683 


— 


— 






— 


14 


21 


56 


106 


-- 


— 






-- 


15 


21 


56 


877 


— 


— 






-- 


IBM 


















1 


7 


7 


11 


7 


8 






18 


2 


7 


7 


30 


10 


17 






21 


3 


3 


4 


63 


2 


13 


48 


1.7 


4 


4 


15 


15 


60 


24 


24 






23 


5 


15 


15 


332 


453 


1078 






154 


6 


31 


31 


341 


60 


417 






— 


7 


12 


50 


391 


F 


74 






• - 


8 


12 


37 


63 


453 


33 








9 


50 


15 


957 


5561 


F 






-- 



(a) Values are from [73. 

(b) Values are from [71. 

(c) Values are from L7J . 

(d) The algorithm presented in this paper. 

(e) The algorithm discussed in [8] , times are from t8] . 

(f) Iterations are the total pivot steps required to 
solve the problem. 

(g) The number of constraints. 

(h) The number of variables, excluding slacks. Each prob- 
lem has as many slacks as it does constraint equations. 

(i) An "F" means that the problem was unsolved after 3000+ 
i terat ions. 

(j) Times to solve problem on IBM 360 series computers. 



Figure 3 
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COMMON A ( 25, 25), B( 25, 25), IB! 25), H( 2 5 ) , AX ( 25) 

COMMON X< 25 ) , XLB ( 25),XUB( 25 ) , X I ( 25), Y( 25), R( 25) 
COMMON XSOL ( 25),IRVl 25) 

COMMON MM,N,M,JS,E,IR,AO,W,Z, IM,DI, IT R , X H, AH , HH, EE , BND 
COMMON ZZ, OJS, ITER,L, IW,NODEO, I NTRE S , XT , XK , I D , RL , BS I GN 
COMMON IREAD, IWRITE, IBNB 
DIMENSION IT I T LE I 1 5 ) 

DATA I STOP / ' ST OP • / 



C PROGRAM TO ENTER DATA AND CALL LP. 

C MM = NUMBER CE CONSTRAINT EQUATIONS. 

C N=NUMBER OF VARIABLES — THE # OF ELEMENTS IN THE X-VECTOR. 

C BND = 0 MEANS NO BOUNDS. 

C BND = 1 MEANS BOUNDS. 

C IL IS A CONTROL VARIABLE TO SIGNAL THE END OF DATA FOR A 
C PARTICLUAR READ STATEMENT. I L= 1 FOR LAST CARD, 0 OTHER. 

C INTRES = C MEANS THAT THE PROBLEM IS A ALL INTEGER 
C RESTRICTED PROBLEM. 

C INTRES = 1 MEANS THAT THE PROBLEM IS A MI X ED- 1 NT EGER 
C PROBLEM. 

C IRV(I) = 0 MEANS THAT THE I-TH VARIABLE IS NOT INTEGER 
C RESTRICTED. 

C IRV(I) = 1 MEANS THAT THE I-TH VARIABLE IS INTEGER 
C RESTRICTED. 

C NOTE THAT THE SL AC K/ SURPLUS VARIABLES DO NOT HAVE TO BE 
C INTEGER RESTRICTED. 

C WHEN I NT RES= 0, THE VALUE OF IRV(I) IS IMMATERIAL, THE 
C PROGRAM WILL HAVE ALL IRV(I)=0, BUT WILL DISREGARD THIS. 

C M IS THE VARIABLE TO AUGMENT THE MARTI X FOR THE OBJ. FUNC. 
C H ( I ) ARE THE RIGHT HAND SIDES, HIM) = INITIAL Z. 

C H(M)=INITIAL Z, WHICH IS ZERO UNLESS OTHERWISE PUT IN PRO. 
C A I M, J ) ARE THE COSTS. 



IWRITE = 8 
IREAD = 4 

2 FORMAT !2I5,F5.0, 15,15A4) 

4 FORMAT I 1 1, 21 3,F5.0 ) 

6 FORMAT ( II ,I3,F5.0,I2) 

10 FORM AT ( 1H ,10F8.0> 

14 FORMAT (1H0»10X, • HAVE FINISHED '/'l') 

49 FORMAT I 1H 1 , 3 5X , 1 5A4// ) 

51 FORMAT I • MINIMIZE: • ) 

52 FORMAT !• SUBJECT TO:') 

53 FORMAT I 1X,24F 5.0) 

54 FORMAT I ' + ' , 12 4X, ' = • , F 5.0 ) 

55 FORMAT (• THIS PROBLEM HAS', 15,' CONSTRAINT EQUATIONS, 
1 ,15,' VARIABLES, AND NO BOUNDS',/) 

56 FORMAT (• THIS PROBLEM HAS', 15,' CONSTRAINT EQUATIONS, 
1 ,15,' VARIABLES, AND HAS UPPER BOUNDS',/) 

57 FORMAT! • WITH THE FOLLOWING UPPER BOUNDS ON THE VARI', 
l'BLES (A -1 INDICATES NO BOUND ON THAT VARIABLE):') 

58 FORMAT!' THE FOLLOWING VARIABLES MARKED WITH A ONE A', 
l'PE INTEGER RESTRICTED:') 

59 FORMAT I IX, 2415) 

15 READ ! IREAD, 2 ) MM , N , BND , I NTRE S , I T I TL E 
IF ( ITITLE(l).EQ.ISTOP) STOP 

M = MM + 1 
DO 1 1=1, M 

hi n=o . 

DO 1 J=1,N 
1 A! I , J I =C . 

DO 9 J=1,N 
I R V { J ) = 0 
XUB! J I =- 1 . 

9 X LB ! J ) =0 . 

3 READ ! IREAD, 4) IL, I ,J,A! I , J) 

I F ( I L. EQ.O )G0 TO 3 

5 READ 1 IREAD, 6) IL, I ,H1 I ) 

I F ! I L . EQ.OIGO TO 5 
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7 READ ( IR EAO, 6 ) I L ,J , A ( M , J ) , I R V ( J ) 
IMIL.EQ.OGO TQ 7 

I F ( BND . EQ . C ) GD TO 8 
11 REA0(IREAD«5 ) IL,J,XUB(J) 

IF ( IL.EQ.CIGQ TO 11 

8 WRITE ( I WR IT E » 49 ) ITITLE 

IF (BND.EQ.O.) WRITE (IWRITE«55) MM, N 
■ IF ( BND. EQ . 1 . ) WRITE (IWRITE, 56) MM » N 
WRITE ( I WR I T E » 5 1 ) 

WRITE ( I WR IT E» 53 ) ( A ( M , J ) , J= 1 , N ) 

IF ( BN 0. FQ.O) GO TO 61 
WRITE! IWRITE, 57) 

WRITE ( I WR IT E t 53 ) ( XUB ( J ) » J = 1 » N ) 

61 IF ( INTRES .EQ.O) GO TO 62 
WRITE (IWPITE,5B) 

WRITE ( I WR I TE » 59 ) ( I RV ( I ) , I = 1 , N ) 

62 WRITE! IWP ITE, 52) 

DO 60 L = 1 » MM 

WRITE ( I WR I TE » 53 ) ( A ( L , J ) » J=1 , N ) 

60 WRITE ( IWR ITE, 54) H( L I 
C EXECUTION TIMED FROM THIS POINT 
CALL LP 

WRITE ( I WR ITE, 14 ) 

C END OF EXECUTION TIMING 

WRITE ( I WR I TE r 49 ) ITITLE 
DO 47 J=1,N 
I F ! X SOL ( J ) )48,47,48 
48 WRITE! IWRITE, 50) J,XSOL!J» 

50 FORMAT (1H0 ,5X,5HXSCL( ,I2,2H)=,F13.4) 
47 CONTINUE 

WRITE ! IWRITE, 45)ZZ 
45 FORMAT (1HC,3HZZ=,F10.5) 

GO TO 15 
END 



SUBROUTINE LP 

COMMON A! 25, 25), B( 25, 25), IB! 25), H( 25), AX! 25) 

COMMON X! 25 ) , XLB ( 25), XUB! 25), XI! 25), Y( 25), R( 25) 
COMMON XSOL ( 25) , I R V ( 25) 

COMMON MM»N,M, JS,E, IR,AD,W»Z» IM,DI , I TR , XH , AH ,HH , EE , BND 
COMMON ZZ»DJS» IT ER , L , IW,NODEO, INTRE S , XT , XK , I D , RL , B S I GN 
COMMON IPEAD, IWRITE, IBNB 



C THIS IS THE SUBPROGRAM WHICH CONTROLS THE SIMPLEX 
C PROCEDURES, THE INITIAL SET-UP OF THE COMPUTATIONAL 
C ARRAYS, ANC THE TERMINATION OF THE PROBLEM. 



NODEO = 0 
ZZ=1000000 . 

IW=C 
Z=H( M) 

L=0 

M=MM+2 
I D=C 

DO 3 J=1 , MM 
DO 2 1=1, M 

2 B ( I , J ) =C . 

B 1 M, J ) =-l . 

3 B ( J , J ) = 1 . 

C XI(J)=+1. MEANS INCREASE. 
C X I ! J ) =-l . MEANS DECREASE. 
C XI ( J) = 2. MEANS BA SIC. 

C XI ( J ) =0 . MEANS NONBASIC. 
DO 24 J= 1 , N 
XI ( J )=0. 

24 X ( J ) =0 . 

DO 4 1=1 , MM 
Y ( I )=H( I ) 

4 IB ( I )=-I 
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AD=1 . 

D I = 1 . 

E= .5 
W=0. 

DO 5 J =1 , N 
A ( M, J ) =0 • 

DC 5 1=1, MM 

5 A(M,J)=A(M,J)-A( I , J) 

DO 6 1=1, MM 

6 W=W+H( I ) 

C W = SUM OF THE H(I>, 1=1,. ...MM. 

WRITE! IWRITE,14)W 

14 FORMAT ( 1H0 »5X , ' INITIAL W= • ,F12.4) 

I TR=C 

I T ER=0 

8 JS=0 

IF ( ITER. GT. 225) WRITE ( I WR I TE , 51 00 ) 

5100 FORMAT ( 1 C X I TERATIONS ARE ABOVE LIMIT SET IN PROG.') 
IF ( ITER. GT. 225) RETURN 
DO 9 J = 1 , N 

I F { X I ( J ) • E Q. 2 . ) GO TO 9 
D= 0 . 

DO 7 K = 1 , MM 

7 D=D+B( M, K) *A( K, J) 

IFIW.EG.O. ) D=D+A ( M , J ) 

IF(ABS(D).LE.E)GO TO 9 
IF (D )15,9, 16 

15 I F ( XUB ( J ) . EQ . - 1 . ) GO TO 19 
IFIABS (XUB(J)-X(J ) ) .LT • E ) GO TO 9 

C INCREASE. 

19 XI (J)=l. 

GO TO 17 

16 IF(ABS( X{ J )-XLB( J ) ) .LT.EIGO TO 9 
C DECREASE. 

XI ( J ) = -l. 

17 IFIJS.GT.C )G0 TO 10 

61 JS=J 

D J S=D# X I ( J ) 

GO TO 9 

10 IF(DJS-D*XI< J ) .GE.E)GO TO 61 

9 CONTINUE 

IF( JS. EQ.C )G0 TO 101 

WRITE! I WRITE, 2 13) JS,DJS,XI< JS) 

213 FORMAT ( 1H , 5X , 3H J S= , I 5 , 5X ,4HD JS= , F 10 . 4 , 5X , 3HX 1 = , F 1 0 .4 ) 
221 CALL ELEM 

I TER= I TER+ 1 
IF (L.LT.C ) GO TO 529 
IF! IR.GT .C )GU TO 36 
IF (BND.EQ.C)GO TO 33 
IF ( BND.EQ.l.AND.L.EQ.O) GO TO 33 
IF(XH. EQ. 1000000. ) GO TO 33 
CALL FLAG (1105) 

36 DO 38 1=1, MM 

I F ( I B ( I ) >39,528,40 

39 Y< I)=Y(I )-AX( I )*XH*XI( JS) 

GO TO 38 

40 K= IB ( I ) 

X(K)=X(K ) — AX ( I )*XH*XI < JS) 

38 CONTINUE 

X(JS)=X(JS)+XH*XI ( JS) 

IF(W.EQ.C. )G0 TO 62 
W=W+AX(M)*XH*XI!JS) 

62 Z=Z+AX(MM+1)*XH*XI (JS) 

WRITE! IWRI TF ,28) W, Z, ITER 

28 FORMAT ( 1 H , 10X , 2HW=, F 10 .4 , 10X , 2HZ= , F 10 .4, 5X , 5H I TER= , I 5 
IF ( IR) 528, 50,44 
44 CALL PIVOT 

CALL FLAG ( 3083 ) 

50 IF! IR.EQ.O) XI ( JS)=0.0 
IF(W.GT.E)GO TO 8 
IF! I W. EQ.l )G0 TO 3 

w=o . 
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WR ITF ( IWRI TE ,73) 

73 FORMAT (1H0,10X, • START PHASE TWO •) 

GO TO 8 

33 WRITE! IWRITE.704) 

WRITE (6,5701) L,R(L),Z,ZZ 
704 FORMAT ( 1HC ,10X ,18H UNBOUNDED SOLUTION) 

RETURN 

101 I F ( W .GT . F ) GO TO 35 
WRITE! IWRI TE ,204) 

204 FORMAT (1HC ,5X,6HRESULT ) 

DC 31 J=1 » N 

I F ( A B S ( X ( J ) ).LT.E)GO TO 31 
WRITE! IWRITE,200) J,X( J) 

200 FORMATdH , 5 X , 2HX ( , I 2 , 2H ) = , F 10 . 4 ) 

31 CONTINUE 

Z=SIGMAINT(DI*ABS(Z) + .5)/DI,Z) 

WRITE! 1WRITE,527)Z 
527 FORMAT (1H0,5X,2HZ=,F17.5) 

C HAVE CONTINUOUS SOLUTION. 

IF ( INTRES.EQ. 1) GO TO 301 
IF (Z+l-E.GE.ZZ) GO TO 134 
301 IF (Z+E.GE.ZZ) GO TO 134 
I F ( D I- 1 . ) 6C1 , 60 1 , 603 
601 WRITE! IWP I T E , 70 1 ) 

WRITE (6,5701) L,R(L),Z,ZZ 

5701 FORMAT ( 5X , • R ( • , I 3 , • ) = • , F 1 4. 6 , 5X , • Z = ',F11.4, 

1 5X , * ZZ = • , F 1 1 .4,/ ) 

701 FORMAT ( 1HC, 5X .22HHAVE INTEGERS REQUIRED) 

GO TO 51 

603 IF ( INTRES ,EQ. 1 ) GO TO 550 

IF ( A B S( ADS! Z)-AINT(ABS(Z)+E) ) .GE.E) GO TO 6032 
550 DO 637 1=1, MM 

IF ( INTRES.EQ. 1. AND. IRV! IB! I) ) .EQ.C) GO TO 637 
IF! IB! I) )637,637, 33 
83 J= I B ( I ) 

IF(ABS(X(J)-AINT (A8S(X(J) )+E) ) -E ) 637, 637, 6032 
637 CONTINUE 

C HAVE ALL INTEGERS REQUIRED. 

WRITE! I WRITE, 701 ) 

WRITE (6,5701) L,R(L) ,Z,ZZ 
51 DO 45 J=1,N 

XSOL(J) = AINT (X(J) + 0.5) 

45 IF ( INTRES.EQ. 1. AND. IRV! J) .EQ.O) XSOL(J) = X(J) 

ZZ = Z 

IF (BND.EQ.O) GO TO 46 

IF ( BND.EQ.l. AND. L. EQ.O) GO TO 46 

GO TO 34 

6032 WRITE! IWPITE.703) 

WRITE (6,573 1)L»R(L)»Z,ZZ 
703 FORMAT!* HAVE FRACTIONS') 

BND= 1 . 

IF! ID. EQ. 500 ) GO TO 41 
42 CALL BOUND 
GO TO 8 

35 WRITE! I WRITE, 205) 

WRITE (6,5701)L,R(L),Z,ZZ 

205 FORMAT! 1H0, • INFEASIBLE*) 

IF (BNC.FQ.O. ) GO TO 528 

GO TO 34 

134 WR ITE ( I WR ITE, 702 ) 

WRITE (6,5701) L,R(L),Z,ZZ 

702 FORMAT (• Z VALUE GREATER THAN SOLUTION UPPER BOUND •) 

34 CALL UNBND 

IF (L.EQ.l .AND.N0DE0.EQ.2 ) GO TO 46 
IF(L.FQ.1.AND.ABS(R(L) ).GT.500) NODEO = 2 
GO TO 8 
K = XK 

IF(XI(K) . EQ. 2. ) GO TO 42 
IF(ABS (XT— X(K ) ) — E) 260,210, 220 
220 I F ( R ( L ) ) 230,8 ,240 
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240 IF (XT.LE.X(K)) GO TO 250 
JS = K 
XI (K) = + 1 . 

L = L - 1 
GO TO 221 
250 ID = C 

XLB(K) = XT 
GO TO 42 
210 I D=0 

XUBI K ) =XT 
GO TO 42 

230 I F ( XT . GE . X ( K ) ) GO TO 210 
JS=K 

XI <K)=-1 . 

L = L - 1 
GO TO 221 

260 IF ( R ( L ) >210,8,250 
529 WRITE (6,5701) L,R(L>,Z,ZZ 
WRITE ( IWR ITE, 5702 ) 

5702 FORMAT ( ' ERROR, L IS NEGATIVE.') 

46 00 47 J=1,N 

I F ( X SOL ( J ) )48, 47,48 

48 WRITF< IWRI TE,49) J,XSCL( J) 

49 FORMAT (1HC ,5X,5HXS0L(,I2,2H)=,F13.4) 

47 CONTINUE 

WRITE ( IWRITE, 52) ZZ 
W R IT E ( 6, 52 I ZZ 

52 FORMAT (1HO,3HZZ=»F10.5) 

528 RETURN 
FND 



SUBROUTINE AXE 

COMMON A ( 25, 25), B( 25, 25 ) , I B ( 25), H( 25>,AX( 25) 

COMMON X( 25 ) , XL B ( 25),XUB( 2 5 ) , X I ( 25), Y( 25), R( 25) 
COMMON XSOL ( 25),IRV( 25) 

COMMON MM,N,M»JS»E, IR,AD»W»Z»IM»D1» I TR , XH , AH , HH, EE ,BND 
COMMON ZZ,OJS,ITER,L, IW,NOOEO, I NTRES , XT , XK , I D, RL , BS IGN 
COMMON IREAD, IWRITE, I8NB 



C THIS IS THE SUBPROGRAM WHICH COMPUTES THE A ( JS ) COLUMN 
C NEEOED TO PERFORM THE PIVOT STEP. 



DO 21 1=1, M 
AX ( I ) = 0 • 

DO 21 K=1,MM 

AX (I )= AX ( I ) + B ( I » K ) *A ( K , JS ) 

21 CONTINUE 

A X ( MM+ 1)=AX(MM + 1 ) + A ( MM+ 1 , J S ) 

RETURN 

END 



SUBPOUT INF PIVOT 

COMMON A ( 25, 25), B( 25, 25 ) , I B( 25), H( 25),AX< 25) 
COMMON X( 25 ) , XLB ( 25),XUB( 25 ) , X I ( 25), Y( 25), R( 25) 
COMMON XSOL ( 25) , IRV< 25) 

COMMON MM, N, M, JS, E, IR,AD,W, Z, IM,DI , ITR.XH, AH , HH, EE, BND 
COMMON ZZ»DJS,ITER»L, I W , NODEO , I NTRES , X T , XK , I D, RL , BS I GN 
COMMON IR EAD, IWR ITE, IBNB 



C THIS IS THE SUBPROGRAM WHICH DOES THE SIMPLEX PIVOT STEP. 



AD=AD*ABS( AX( IR ) ) 

DI=AINT(DI*ABS(AX( IR) ) + .5) 

84 WRITE( IWRITE, 214) IR,AX(IR) 

214 FORMAT ( 1H0 » 2 5X , 3H I R= , I 3 , 3X , 3HAX= , F8. 4 ) 
IF(IB(IR))1,1,2 
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2 K = I B I I R ) 

XI (K )=C. 

I I e ( I R ) = J S 
XI I JS>=2. 

E=.5/DI 

DO 29 i = i,m 

IF ( I .EQ.IR)GO TO 29 

P = -AX ( I ) / AX ( I R ) 

H ( I ) =H ( I ) + P*H I IR) 

DO 31 J=1,MM 

31 B( I, J)=B( I ,J )+P*B< IR» J ) 

29 CONTINUE 

P= 1 . / AXI I R ) 

H ( IR ) = H( IR ) *P 
DC 33 J=1,MM 
33 B ( IR , J ) =B I IR,J)*P 
91 DO 92 1=1, M 

HI I )=SIGN< A I NT IDI*ABS(H( I ) )+.5)/DI ,H{ I ) ) 

DO 92 J= 1 , MM 

9 2 PI l, J)=S IGMAINTI DINARS (B( I,J) )*.5)/DI,B< I,J) ) 

W = AINTIDI*W+0.5)/DI 
Z = SIGN(AINT(Dl*ABS(Z) + 0.5)/DI,Z) 

95 I T R= I T R+ 1 

WRITE! I WRITE, 2 10) I TR , W , Z , DI , AD 
210 FORMAT I 1HC ,5X ,4HITR=, I 5 , 5X , 2HW= , F 1 2. 8 , 5X , 2HZ= , F 1 2 . 5 , 5X 
1 ,3X,3HAO=, FI 2 . A ) 

RETURN 

END 



SUBROUTINE VAR 

COMMON A I 25, 25), B( 25, 2 5 ) , I B I 25), HI 25),AX( 25) 
COMMON XI 25 ) , XL B I 25),XUB< 25), XII 25), Y( 25), R( 25) 
COMMON XSOLI 25),IRVI 25) 

COMMON MM, N,M, JS, E, IR , AD, W, Z, IM, DI , I TR , XH , AH , HH, EE , BND 
COMMON Z Z, DJS , ITER , L , IW,NODEO, I NTRES , XT , XK , I D, RL , BS IGN 
COMMON IREAD, IWR ITE, IBNB 



C THIS SUBROUTINE PICKS OUT THE NON-INTEGER VARIABLE IN THE 
C CURRENT SOLUTION THAT WILL BE BOUNDED AND BRANCHED UPON. 



XLDU=0 .G 
XLDL=0.0 
I ML=0 
I ML=C 

62 DO 5C IV = 1 , MM 
XDU= 1CC0C0C 
XDL= 1000000 

ir ( I B I IV) ) 50,50,55 

55 IF I IRVI I B I IV) ) .EQ.O.AND. INTRES.EQ. 1 ) GO TO 50 
IFIABSIXI I B( IV ) )-AINT( XI IBI IV) )+E) )-E ) 50, 50, 60 

60 DO 61 JV=1 , MM 

WRITE I I WR I TE » 44 ) B I M , J V) , B I l V , JV ) 

44 FORMAT (• BI M , JV ) = • , F 1 4 . 4, • , B I I V , J V ) = • , F 9 . 4) 
IFIABSIBI M,JV) ) .LT.E) GO TO 61 
IFIABSIBI IV, JV) ) .LT.E) GO TO 61 
I F I B I I V, JV ) ) 64,61,63 

63 XDUJV=ABSt BIM, JV) /BI IV, JV) ) 

XDL JV= O.C 

I FIXDUJV.GE.XDU) GO TO 61 
XDU=XDUJV 

WRITE I IWR ITE, 41 ) XDLJ V, XDUJ V , JV, I V,X< IBIIV)),IB(IV) 
GC TO 61 

64 XDLJV=ABS(BIM,JV)/B!IV,JV)) 

XDUJ V= 0.0 

IFIXDL JV.GE.XDL) GO TO 61 
XDL=XDL J V 

WRITE I IWR ITE, 41 ) XDLJ V, XDUJV, JV, I V,X( IB! IV)) ,IB( IV) 

61 CONTINUE 

52 FIDLO = XIIBIIV))- A I NT I X I I B I I V) ) ) 
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53 XLDLIV=XDL*F IDLO 

XLDUIV=XDU*( l.O-F IDLO) 

IF (XLDLIV. IE .XLDL ) GO TO 66 
XLDL=X LDL I V 
I ML= I V 

66 IF(XLDUIV.LF.XLDU) GO TO 50 
X L DU=X LDU I V 

I MU= I V 
50 CONTINUE 

IF (XLDU-XLDL > 68,67,67 

67 I M= I MU 

BS IGN= 1.0 
GO TO 45 

1T0 FORMAT (FI. 01 

68 I M= I ML 
BSIGN=-1.0 

45 WRITE! I WRITE, 10) IM 
WRITE (6 , 10 ) I M 

READ ( 5 , IOC ) DS TOP 
IF (DSTOP.EQ.l .0) STOP 
10 FORMAT ( 1HC , 10X , • FROM VAR I M= ' ,12) 

WRITE (6 ,42) BSIGN 

WRITE ( I WR I TE , 40 ) F I OLO , XDU , XOL , XL DU , XLDL , I MU , 

1 I ML , IM 

40 FORMAT ( • F 1 0L0= • , F9 . 4 , • , XDU= • F9 . 4 , • , X DL= • , F9 . 4 , 

114, ' ,XLDU= ' » F9 .4* , ,XLDL=',F9.4,',I MU= ' ,14, 

2' ,IML=',I4,' , I M= ' , 14,/) 

41 FORMAT ( ' XDL JV=' ,F9.4, • ,XDUJV=' ,F9.4, • ,JV=' , 14, • , IV=' , 
1I4,',X(IB(IV))='»F5.1,'»IB(IV)='»I4»/) 

42 FORMAT!' B S I GN= • , F9 . 4 ) 

IF ! IM.EO.O )STOP 
RFTUPN 

END 



SUBROUTINE REDO 

COMMON A! 25. 25), B( 25, 25), IB! 25), H! 25), AX! 25) 
COMMON X! 25 ) » XLB ( 25),XUB( 25), XI! 25), Y! 25), R! 25) 
COMMON XSCL! 25),IRV( 25) 

COMMON MM,N,M, JS,E, IR, AD,W,Z, IM,OI , ITR,XH, AH,HH, EE , BND 
COMMON ZZ,DJS, ITER,L, I W , NODFO , I NTR E S , XT , XK , I D , RL , B S I GN 
COMMON IREAD, IWRITE, IBNB 



C THIS SUBROUTINE REMOVES THE LAST BOUND SO THAT BACK TRACK- 
C ING CAN START, AND RESETS ALL OF THE OTHER BOUNDS. 



5 WRITE (IWRITE, 14) 

14 FORMAT (1H0»10X, • FROM REDO ') 

DO 21 K= 1 , N 
XUB ( K ) =-l . 

XLB! K ) =0 . 

21 CONTINUE 
LM=L-1 

IF (LM.EQ.O) GO TO 17 
DO 18 JT=1,LM 
CALL UNPAK (JT) 

K= XK 

IF ( R ( JT ) ) 19,18,20 

19 XLB! K ) =XT 

WRITE! IWRITE, 11)K,XLB(K) 

11 FORM AT ( 1H0 ,10X,2HK=, I 2 , 5X , 4HXL 8= , F 10 . 4 ) 
GO TO 18 

20 XUB! K ) =XT 

WRITE! IWRITE,12)K»XUB(K) 

12 FORMAT ( 1H , 10X , 2HK = , I 2, 5X, 4HXUB=,F 10. 4 ) 
18 CONTINUE 

17 JT = L 

CALL UNPAK (JT) 

I D=500 

IF(BSIGN.EQ.-1. ) XT = XT— 1 
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IFlBSIGN.EQ.+l.) XT= XT 1 
CALL PACK 
BSIGN = - BSIGN 
R(L) = -R(L> 

WRITE(IWRITE,15)L,R(L) 

15 FORMAT 11HC ,10X,2HL=,I2,5X,2HR=,F14. 8, • XXXX XX XXXXXXX » ) 
1 RETURN 
END 



SUBROUTINE PACK 

COMMON A ( 25, 25), B( 25, 25),IB( 25), H( 25),AX( 25) 
COMMON X( 25 ) » XLB ( 25),XUB( 25 I , X I C 25), Y( 25), R( 25) 
COMMON XSOL ( 25>,IRV( 25) 

COMMON MM, N, M, JS,E ,IR, AO,W,Z, IM,OI , ITR,XH, AH,HH, EE , BNO 
COMMON ZZ, DJS, ITER,L, IW,NODEO, INTRES,XT, XK,I D,RL, BSIGN 
COMMON IREAD, IWRITE, IBNB 



C THIS SUBROUTINE COMBINES INTO A SINGLE NUMBER THE VARIABLE 
C BOUNDED, THE VALUE OF THE BOUND, AND WHETHER IT IS AN 
C UPPER OR LOWER BOUND. 



IF ( XT . EQ. 0 . ) GO TO 1 
I L= 1 
I T=XT 

30 IF ( IT/ 1C) 11,11,20 

20 I L= I L+ 1 
I T= I T/ 10 
GO TO 30 
11 TL= I L 

R ( L ) = ( XK+TL/1C •♦XT/10»**(TL + 1. ) +FLOAT ( 10) ) *! BSIGN) 
WRITE! IWRITE,10)XT,IL,TL,L,R(L) 

10 FORMAT ( 1 HO , 1 C X , • FROM PACK • , 5X , 3HXT= , F 1 C . 4 , 5X , 3HI L= 
1,F IG.4,5X,2HL=,I5,5X,2HR=,F14.8) 

GO TO 2 

1 R ( L) =( XK + FLOAT ( ID) )*( BSIGN) 

2 RETURN 
END 



SUBROUTINE UNPAK (JT) 

COMMON A! 25, 25), B( 25, 25), IB! 25), H( 25), AX! 25) 

COMMON X! 25), XLB! 25),XUB! 25), XI! 25), Y( 25), R! 25) 
COMMON XSOL! 25),IRV< 25) 

COMMON MM,N,M, JS, E, IR,AD,W,Z, IM,DI , I TR , XH , AH , HH , EE , BND 
COMMON ZZ, DJS, ITER, L, IW,NODEO, INTRES , XT , XK , I D, RL , BS IGN 
COMMON IREAD, IWRITE, IBNB 



C THIS SUBROUTINE RETREVIES THE INFORMATION COMBINED IN 
C SUBROUTINE PACK. 



RL=ABS(R! JT) ) 

I D=0 

IF {RL.LT.5C0) GO TO 10 
RL=RL- 500. 

I D=500 
10 IH=RL 
XK=IH 

T L= ( RL-XK ) *10 . 

I L=TL 
RX= I L 

IXT= (TL-RX )*10 .* : <'RX+. 5 
XT=I XT 

BSIGN = SIGN!1.0,R(L) ) 

RETURN 

END 
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SUBROUTINE BOUND 

COMMON A ( 25, 25), B( 25, 25), IB! 25), H( 25),AX( 25) 

COMMON X( 25 ) , XL B t 25),XUB( 25), XI! 25), Y( 25), R( 25) 
COMMON XSOL ( 25),IRV( 25) 

COMMON MM,N»M, JS,E , IR,AD,W,Z, IM,DI , I TR.XH , AH , HH , EE , BND 
COMMON ZZ, DJS, ITER,L, IW,NODEO, INTRE S , XT , XK , I D , RL , B S I GN 
COMMON IREAD, IWRITE, IBNB 



C THIS SUBROUTINE SETS THE VALUE OF BOUND OF THE NEXT BRANCH 



IDBND = ID 
LSET = C 

I F ( IO.EQ.500) GO TO 2 
L = L+1 
LSET = 1 
23 CALL VAR 

IE ( IM.GT.O) GO TO L7 
WRITE (IWRITE, 14) IM 
L = L - 1 
RETURN 
17 K= I B ( I M ) 

I h=X (K ) *D I +. 5 
XT = I H 
X( K)=XT/DI 
I H=X ( K )+E 
32 XT= I H+ 1 

IF(BSIGN.EQ.-1. ) GO TO 31 
CALL PACK 
GO TO 6 

31 Y ( I M ) = XT— X ( K ) 

X ( K ) =XT 
XLB( K) =XT 
I B ( I M ) =- I M 
M=M+ 1 
XI ( K ) =0. 

DO 1 J =1 , MM 
B(M, J)=B( IM, J) 

1 B ( IM, J)=-B< IM, J) 

H ( I M ) =— H ( I M) 

W=Y ( IM ) 

I W=0 

WRITE! IWRITE, 16) W 

I F ( I DBND.EQ.500 ) RE TURN 

XK=K 

CALL PACK 
GO TO 3 

2 K= XK 

IF(ABS(XT— X(K) ) — E ) 21,21,22 

21 I D=C 

I F ( B S I GN.LT.O.O) GG TO 30 
XUB( K ) =XT 
GO TO 23 
30 XLB(K)=XT 
GO TO 23 

22 IF(XT.GE.X (K) . AND . BS I GN . EQ . + 1 . ) GO TO 21 
IF(XT.LE.XtK) .AND.BSIGN.EQ.-l. ) GO TO 21 
DO 5 1=1, M 

I F ( IB( I) . E Q. K ) GO TO 4 
GO TO 5 

4 IM = I 

WR ITE ( IWRITE, 14) IM 

14 FORMAT ( 1H0 , 10X, • FROM BOUND I M= • ,15) 

GO TO 8 

5 CONTINUE 
WRITE! IWRI TE .15) 

15 FORMAT (1H0 ,10X , • ERROR • ) 
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8 ID = 0 

IF (BSIGN.LT. 0.0) GO TO 31 

6 Y ( I M ) = X ( K ) -XT 
X( K) =XT 

XUB ( K ) =XT 

WR I T E ( IWRITE ,10)K,X!K) 

10 FORMAT ( 1H0 ,10X , • FROM BOUND ' , 5X ,2HK= , 12 ,5X ,2HX=, F10 

I B ( I M ) =— IM 
M=M+1 
X I ( K 1 = 0. 

DO 7 J = 1 , MM 

7 B( M, J) =-B( IM, J ) 

W=Y( IM) 

IW=0 

WRITE! I WR I TE » 1 6 ) W 

16 FORMAT ( lHOtlOX, • FROM BOUND W= * ,F10.8) 

IF (LSET.FO.l )RETURN 
L = L + 1 
XK = K 
CALL PACK 
3 RETURN 
END 



SUBROUTINE UNBND 

COMMON A ( 25, 25), B( 25, 251, I B ( 25), H( 25),AX( 25) 
COMMON X( 25 ) , XL B ( 25),XUB( 25), XI( 25), Y( 25), R( 25) 
COMMON XSOL ( 25),IRV( 25) 

COMMON MM , N, M, JS , E , I R , AD , W, Z , I M , D I , I TR ,XH , AH , HH , E E , BND 
COMMON ZZ, DJS, ITER,L, IW,NODEO, I NTRE S , X T , XK , I D , RL , BS I GN 
COMMON IPEAD, IWRITE, IBNB 



C THIS SUBROUTINE RESETS ALL THE BOUNDS AFTER EACH BACKTRACK 
C AND DETERMINES WHICH NODE TO BACKTRACK TO. 



6 IF(ABS(R(L))-500.0) 1,1,2 

1 CALL REDO 
I D=5C0 
RETURN 

2 IF(L.EQ.l) CALL REDO 
IF(L.EC.l) RETURN 

L = L - 1 
GO TO 6 
ENC 



SUBROUTINE ELEM 

COMMON A ( 25, 25), B( 25, 25 ) , I B( 25), H( 25),AX( 25) 
COMMON X ( 25 ) , XLB ( 25),XUB( 25), XK 25), Y( 25), R( 25) 
COMMON XSOL ( 25),IRV( 25) 

COMMON MM, N,M, JS, E, IR,AD,W, Z, IM,DI, ITR,XH, AH,HH, EE, BND 
COMMON ZZ,DJS, ITER.L, IW,NODEO, I NTRES , XT , XK , ID,RL,BSIGN 
COMMON IREAD, IWR ITE, IBNB 



C THIS IS THE SUBPROGRAM WHICH OBTAINS THE PIVOT ELEMENT 
C FOR THE GIVEN COLUMN DETERMINED EITHER IN *LP' OR 
C IN 'BOUND*. 

C IT GIVES WHICH VARIABLE WILL LEAVE THE BASIS. 



I R=C 

CALL AXE 
EE=E 

WRITE! IWRITE, 43)E 

4 3 FORMAT ( 1H0 , 10X , 12HFR0M ELEM E=,F10.4) 
IF (BND.EQ. 0. )G0 TO 12 
IF ( X I ( JS).EQ.1.)G0 TO 1 
GO TO 2 

1 XH = XUB(JS)-X( JS) 
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IF (XUB(JS) .EQ.-l. )XH= 1000000. 
GO TO 12 

2 XH=X( JS)-XLB( JS) 

12 DO 25 1=1, MM 

IF (BNC.GT.C. ) GO TO 810 
IF ( AX { I ) .LE.EIGO TO 25 
AH=AX ( I ) 

Hh=H ( I I 

IF(IR)25,811,812 
311 I R = I 

XH=HH/AH 
EE=E/AH 
GO TO 25 
812 CALL MIN (II 
GO TO 25 

810 IF(ABS(AX( I) I.LT.EIGO TO 25 
IF (XI ( JS) . EQ. 1 . IGO TO 8 
C DECREASE. 

I F ( I B ( I) 13,25,4 

4 K= IB ( I ) 

I F ( AX ( I ) 15,25,6 
6 IF (XUB(K) . FQ.-l. IGO TO 25 
HH = XUB ( K I — X ( K) 

AH = — AX ( I ) *XI ( J S ) 

CALL MIN (I I 
GO TO 25 

5 HH=X (K I-XLB(K) 

AH= AX ( I)*X I( JS I 
CALL MIN (I) 

GO TO 25 

3 IF ( AX ( I I 115,25 ,16 
16 HH=Y ( I I 

AH=— AX ( I)*XI(JS> 

CALL MIN (I I 
GO TO 25 
1 5 HH=Y ( I I 

AH=AX( 1 1 *X I ( JS I 
CALL M IN (I I 
GC TO 25 
C INCREASE. 

8 IF ( I B ( I 1)30,25,40 
40 K= IB ( I I 

IF ( AX( 11)6,25,5 
30 I F < A X ( I ) ) 1 6, 25 , 1 5 
25 CONTINUE 
RETURN 
END 



SUBROUTINE MIN ( I ) 

COMMON A ( 25. 25). B( 25, 25 ) , I B ( 25), H( 25),AX( 25) 

COMMON X( 25 ) , XL B ( 25),XUB( 25),XI( 25), Y( 25), R( 25) 
COMMON XS CL ( 25 ) , I RV ( 25) 

COMMON MM,N,M»JS,E»IR,AD,W»Z» IM,OI , I TR , XH, AH ,HH, EE , BND 
COMMON ZZ, CJS, ITER,L, IW,NODEO, INTRES, XT, XK , I D,RL ,BSIGN 
COMMON IREAO, I WRITE, IBNB 



C THIS IS THE SUBROUTINE USED TO CALCULATE THE MINIMUM THETA 



IF (ABS (HH-AH*XH)-EE 170,70,2 

2 IF (HH-AH+XH) 24,70,25 
24 XH=HH/AH 

3 I R = I 
EE=E/AH 
GC TO 25 

7 C IF ( IR. EQ.O ) GO TO 25 
I F ( I B( I ) )74, 25 ,73 

73 IFUR.EQ.O ) GO TO 25 
IF ( IB ( IR) ) 25,25,4 

74 IFUR.EQ.O) GO TO 3 
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IF ( IBdRI ) 4,25,3 
4 IF (ABS( AX I IR) )-ABS(AX( 1)1)25,25,3 
25 RFTURN 
END 
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ASSEMBLY LANGUAGE TIMING ROUTINE 



SETIME START 

SAVE 114,12) 

LR 12,15 

USING SETIME,12 
ST 13, TEMP 

LA 13, SAVE 

STIMER TASK, FIXUP, TUINTVL=TIME 
L 13, TEMP 

RETURN (14,12) »T 
DROP 12 

FIXUP SAVE (14,12) 

RETURN (14,12) 

ENTRY WATIME 

hAT I ME SAVE (14,12) 

LR 12,15 

USING WATIME,12 
L 2,0(1) 





ST 


13 , TEMP 




LA 


13, SAVE 




TTIMER 


CANCEL 




L 


1 , TIME 




SR 


1 ,0 




ST 


1,0(2 ) 




L 


13, TEMP 




RETURN 


( 14 ,12) ,T 


SAVE 


OS 


18F 


TEMP 


DS 


F 


TIME 


DC 


X '08000000 




ENC 
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